Structure of surface vortex sheet between two rotating 3 He superfiuids 
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We study a two-phase sample of superfluid 3 He where vorticity exists in one phase ( 3 He-A) but 
cannot penetrate across the interfacial boundary to a second coherent phase ( 3 He-B). We calculate 
the bending of the vorticity into a surface vortex sheet on the interface and solve the internal 
structure of this new type of vortex sheet. The compression of the vorticity from three to two 
dimensions enforces a structure which is made up of i-quantum units, independently of the structure 
of the source vorticity in the bulk. These results are consistent with our NMR measurements. 

PACS numbers: 67.57.Fg, 47.32. Cc, 05.70.Fh 



Consider an interface separating two phases whose or- 
der parameters are coherent across the boundary. How 
do topological line or planar defects behave when they 
meet the boundary? The coherence rules out simple ter- 
mination at the interface. The remaining alternatives are 
that the defect crosses the boundary or is deflected to 
continue along the boundary. If deflected, how does the 
defect bend onto the interface and what is its structure 
when it lies on the boundary? Such questions are cer- 
tainly relevant for dislocations at coherent grain bound- 
aries in crystals. These questions also appear in liquid 
crystals and in the cosmos 0, • Here we provide an 
answer in the context of superfluid 3 He, where detailed 
results can be achieved by combining experimental and 
theoretical analysis 0,0- 

The important defects in superfiuids are vortex lines or 
sheets, which can be created in a controlled way by rota- 
tion. A crucial property of 3 He superfiuids are multiple 
length scales: The core diameter of a typical vortex is 10 3 
times larger in the A phase than in the B phase. Cor- 
respondingly the vortex energy in the A phase is lower 
and vortices do not easily penetrate from the A to the 
B phase, but form a surface sheet on the phase bound- 
ary (Fig.^). Here we calculate how the vorticity in bulk 
3 Hc-A bends to form such a surface sheet. The calcu- 
lated internal structure of this surface sheet turns out to 
be quite different from the vortex sheet that appears in 
bulk 3 He-A. Finally, we report on NMR measurements 
and show that our calculations are consistent with these. 

To obtain a two-phase sample of superfluid 3 He, the A- 
B interface is stabilized in a gradient of magnetic field H . 
Vortices are created by rotating the sample at angular ve- 
locity Q around the axis z perpendicular to the interface. 
In the A-phase section vortices are formed at a low crit- 
ical velocity w c a so that the average superfluid velocity 
(v s a) approximates solid-body rotation, (v s a) « x r. 
Depending on preparation, the vorticity in the bulk is 
in the form of vortex lines or vortex sheets la, Ml . Both 



structures have a large "soft vortex core" region for which 
the length scale £d ~ 10 /im is set by weak dipole-dipole 
forces. The vortex line is doubly quantized in units of the 
circulation quantum kq — h/2m, where m is the mass of 
a 3 He atom. The bulk vortex sheet has periodic units, 
which also consist of two quanta. The B phase vortices, 
in turn, are singly quantized and have a narrow "hard 
vortex core" with a radius comparable to the superfluid 
coherence length £ ~ 10 nm. Their smaller core radius 
causes the B-phase critical velocity to be at least an or- 
der of magnitude larger. Thus the B phase remains in 
metastable vortex-free state, where its superfluid frac- 
tion is stationary in the laboratory frame: u s b = 0. To 
sustain the difference in superflow velocities at the A- 
B interface, the A-phase vorticity has to curve onto the 
interface where it forms a surface vortex sheet (Fig.^). 

Bending of vorticity into surface sheet: — We calculate 
the macroscopic configuration for bending the vorticity, 
both in the form of vortex lines and bulk sheets. For 
vortex lines we use the Bekarevich-Khalatnikov model . 
The coarse-grained vorticity lj = i V x (v s ) is determined 
from the energy functional 



F 



d 3 j 



<y\u>\ + ((v s ) - v n y 



(1) 



Here p s is the density of the superfluid fraction. The 
first term is the energy of vortex lines, where 7 = 
(2k /tt) In (r v /r c ) (Fig. is on the order of the circu- 
lation of a line. The second term is the energy penalty 
associated with the difference between the average su- 
perfluid and normal velocities, (v s ) and v n — fi x r. For 
vortex sheets we use the functional 
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The latter term is the surface energy of the bulk sheet, 
with the surface tension a. In contrast to Eq. 0), the 
velocity v s = (h/2m)'V(p is calculated exactly (V 2 = 0), 
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FIG. 1: A-phase vorticity (drawn in cross section) curves into 
a surface vortex sheet on the A-B interface (SI ~ 1 rad/s). 



treating the bending sheet as a tangential discontinuity 
of v s . 

In rotation with $1 = Slz, we assume the vorticity to 
occupy the half space z > and not to penetrate to z < 0. 
We select a location far from the axis of rotation so that 
the cylindrical coordinates can locally be approximated 
with cartesian coordinates, taking x in the radial direc- 
tion. We assume that the vortex sheets at z — > oo are 
perpendicular to x You then find that between two 
sheets v = v s — v n = 2Sl(x — c)y, where c is a constant. 
This allows us to write Eq. © as 
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Here b = (3a/ Psfl 2 ) 1 ^ 3 is the separation of two sheets in 
the bulk (z — > oo) and £ their radial displacement. 

Surprisingly we find that both models [JU and (0)] give 
exactly the same form of bending: In spite of different 
physics and approximations, the radial deviation £(z) of 
both vortex lines and vortex sheets is given by 
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For vortex lines a = yj 7/ (40) is on the order of the line 
spacing 2r v . For vortex sheets a = b/y/6 is 0.41 times 
their equilibrium spacing in the bulk. The bending con- 
tour plotted in Fig. becomes horizontal at the radial 
deviation £ = a. After this point the bulk vorticity is 
transformed to a new state of vortex matter - a surface 
sheet. The vorticity in the sheet grows linearly with dis- 
tance r from the center and finally escapes to the vertical 
sample boundary. 

Structure of surface vortex sheet: — The order param- 
eter of superfluid 3 He, a 3 x 3 matrix A^j, takes in the 
A phase the form A^j — Axd^(rhj + ihj). Here rh, h, 
and I are three orthogonal unit vectors that specify the 
orbital part of the Cooper pair wave function. The vec- 
tor d specifies the orientation of the spin part. Here we 



consider only continuous vortex structures where the am- 
plitude Aa is constant, and the circulation arises solely 
from a smooth orientational winding of the triad (rh, h, 
I) in space. The vortex textures are determined from the 
energy functional F — J d 3 r /a + J B d 3 r /b, where 
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v A = v s a — v n and v s a — (h/2m)^ i rhiWhi. The first 
nine terms give the gradient energy, while the last two 
terms are the dipole-dipole and external field energies. 
This functional replaces Eqs. and (01 when the reso- 
lution is increased from r v to the dipolar coherence length 
£d = (h/2m)yj p\\/ \d. At the A-B interface one has 
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as boundary conditions, where s is the normal of the in- 
terface. The rotation matrix R and the phase angle <f> 
are quantities appearing in the B-phase order parameter 
A^j — Ase^Rnj- On the A-phase side the role of the 
phase angle is played by the rotation angle of rh and h 
around I. The boundary conditions © imply the coher- 
ence of the phase angle across the interface. 

We calculate the order parameter in the surface sheet 
by minimizing the total energy F numerically. The main 
assumption is that the solutions are homogeneous in the 
(radial) x direction and periodic in the perpendicular 
(azimuthal) y direction. Locally these assumptions are 
approximately satisfied everywhere except near points 
where new vorticity enters the sheet. The velocities far 
above and below the sheet satisfy — Va\ = 2kq/L v 
(or L y = 2kq/ '(Or)) when there are two circulation quanta 
per period L y . Since the velocity is effectively screened 
on the A-phase side by vorticity, we can take Va ~ 0. 
The experiment is performed in a magnetic field H \\ z 
which locks d _L H ^lj . Since the variation of d has only 
a minor effect, we approximate the B phase R with a con- 
stant. In contrast, the phase field (j>(y, z) on the B phase 
side has to be properly included, via /b = \p s {v S B — v n ) 2 
and u S B = (fi/2m) V0. Otherwise a surface vortex sheet 
is not stable, as is the case for a sheet which would coat 
a solid container wall 0. 

Depending on the density of vorticity in the surface 
sheet, we obtain two different textures which both incor- 
porate ^-quantum vortex units. These textures are in- 
dependent of the initial ansatz, ie. whether a 2-quantum 
vortex line or one period of the bulk vortex-sheet is placed 
above the A-B interface at the start of the iterative en- 
ergy minimization. This means that the bulk structures 



FIG. 2: Orbital texture of the surface sheet at low density 
of vorticity (L y = 8£d) in a frame where v n = 0. The cones 
point along I and their yellow stripes indicate the rotation of 
m and n around I. On the A-B plane (x — y plane) I is parallel 
to the interface (Eq. (JSJ). Its orientations there are shown in 
the diagram on the bottom. The four highlighted regions 
(i«±i) are the centers of ^-quantum cores which pairwise 
form 1-quantum composites. Two ^-quantum cores in one 
pair are separated by d c = 0.26 £d + 0.135 L y - 0.0027 L 2 y /£ d , 
when 5.5 Cd < L y < 20 f d - 

loose their identity and transform at the A-B interface to 
surface sheet textures. 

The low-density texture in Fig. |2 has all the vortic- 
ity aligned on the A-B interface. It separates into two 
composite 1-quantum vortices. These consist of two 
^-quantum cores, although the vorticity V x v$\ (not 
shown) cannot be divided into distinct i-quantum units. 
At high density a different packing of the vorticity be- 
comes energetically favorable. The resulting texture in 
Fig-Elhas two ^-quantum cores on the A-B boundary and 
the remaining circulation localized as a 1-quantum vor- 
tex above the A-B interface. Here the vorticity V x v s a 
is maximized in the ^-quantum regions. The distinguish- 
ing feature of the different vortex components is the solid 
angle which their orientational distribution of I covers. 
For instance, the circular 1-quantum vortex above the 
A-B plane in Fig. [21 includes all orientations of the pos- 
itive hemisphere with respect to the x axis, while the 
two ^-quantum vortices each cover one quadrant on the 
negative hemisphere. 

The first-order transition between the two textures 
takes place at L y = 5.7 or when the velocity difference 
in shear flow — Vj£\ « 2.8 mm/s. This value is in the 
middle of our measuring range. In the calculations the 
transition is hysteretic, especially on moving from high 
to low density. Transitions between the two textures are 
thus expected as a function of Q and r: The density of 
the surface vorticity increases with r as ^Or / ' kq. Thus at 
high f2 one expects to find the low-density texture in the 
center and the high-density one outside a critical radius. 



FIG. 3: Orbital texture of the surface sheet at high density 
of vorticity (L y — 5£d)- With increasing density, the two 
^-quantum cores in the center of Fig. [2] form a 1-quantum 
composite, a circular meron, which pops above the A-B plane, 
as seen in this figure. The highlighted regions on the A-B 
plane, the two outermost ^-quantum vortices from Fig. [5] are 
now isolated, but weakly bound by the meron above the A- 
B interface. The separation between the ^-quantum cores is 
d c = 0.0743 £d+ 0.507 L H , when 2£ d < L v < 7.5£d- 



Transfer of vorticity across A-B interface: — Our NMR 
measurements of the two-phase sample show that the A- 
B interface is stable up to high rotation. Ultimately at a 
critical angular velocity Q c , which corresponds to a crit- 
ical B-phase superflow velocity v c rs fl c R with respect 
to the container wall at r « R, the A-B interface un- 
dergoes an instability and a small number of circulation 
quanta manage to cross the interface and form the first 
vortex lines in the B phase. An analysis of this event as 
a function of temperature and barrier field shows that v c 
is determined by the hydrodynamic rigidity of the inter- 
face, on which the texture of the surface vortex sheet has 
little effect. The upper limit for the density of vorticity in 
the surface sheet is thus placed by the shear-flow instabil- 
ity. With our solenoidal barrier magnet this means that 
£l c < 1.6rad/s. If the rotation is increased above fl c , 
then the instability occurs repetitively at the constant 
critical velocity v c , as analyzed in Fig.0J 

From Fig. 0| it is concluded that the number of circula- 
tion quanta AN, which is transferred across the A-B in- 
terface in one event, can be either odd or even. Measure- 
ments accumulated in the constant conditions of Fig. 01 
on more than 200 instability events display a smooth 
probability density distribution P(AN), which is cen- 
tered around (AN) w 8 and where both odd and even 
values of AA^ are equally likely. This indicates that the 
vorticity from the bulk A-phase breaks through the A- 
B interface as single quanta, although the source object, 
the A-phase vortex, is doubly quantized. Thus after an 
instability event the number of quanta in the B-phase 
section of the sample can be odd, while it is in the bulk 
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FIG. 4: NMR measurements on the shear-flow instability of 
the A-B interface provide information on the A-phase surface 
vortex sheet. (A) NMR absorption spectra of 3 He-B in the 
vortex- free state (N = 0) and with 69 vortex lines (N = 69). 
The ratio of the Larmor and counterflow peak heights, mea- 
sured at constant conditions (at = f2 rc f < f2 c ), is a linear 
function of N. (B) A repetitive sequence of instability events 
with increasing Q. The data points plot the peak height ratio 
from NMR spectra measured at f2 fe f- The spacing of the hori- 
zontal grid lines is the calibrated equivalent of one vortex and 
yields the number of new B-phase vortex lines AN created 
in each instability event (explicitly given at each discontin- 
uous step). The dashed vertical lines (and arrows in panel 
D) denote the il values where the Larmor peak height rises 
abruptly while fi is slowly increased by a small increment. 
The sloping dashed line is a fit through the corner points 
and defines the effective radial location i? c ff of the instabil- 
ity: N = (2w/k )RIs (fl - fic). (C) One of the instability 
events in greater detail. The resolution allows us to conclude 
that AiV is a small random number, which can be even or 
odd. (D) Discontinuous superflow velocity with increasing Q 
at r — R c g: v = \v s b — v n \ = QR c g — KoN/(2wReg). The 
horizontal dashed line is equivalent to the sloping dashed line 
in panel B and defines the mean critical velocity v c - 

A phase even. This means that one unpaired quantum 
must be accommodated and stored in the surface sheet 
in stable state. These properties remain unchanged if the 
A-phase section of the sample is arranged to support the 
bulk vortex sheet, in any of its different global configura- 
tions, as explained in Ref. 0. The measurements support 
the conclusion from the texture calculations that, com- 
pared to the bulk A-phase vortex textures, the surface 
sheet is made up from smaller building blocks and has 
autonomous structure, independently of the bulk vortex 



textures above the A-B interface. The experiments can- 
not as yet distinguish between the two calculated tex- 
tures, but in both cases the vorticity can be combined 
to 1-quantum units which allows the transfer of vorticity 
from the A to the B phase. 

Discussion: — Stable vortex sheets are discussed in co- 
herent quantum systems with a multi-component order 
parameter. So far two examples have been experimen- 
tally identified, both in anisotropic superfluid 3 He-A, 
namely the bulk vortex sheet and the present surface 
sheet. The stability of their structures is based on differ- 
ent principles. The bulk sheet is topologically stable: Its 
linear quantized vorticity consists of 1-quantum building 
blocks, known as merons or Mermin-Ho vortices, which 
come pairwise as a combination of a circular and a hyper- 
bolic unit. These are confined within a domain- wall- like 
planar defect. In the surface sheet the quantized vortic- 
ity is confined to two dimensions by the hydrodynamic 
binding of the circulation to the interfacial boundary. 
The texture is made up of ^-quantum building blocks 
to satisfy the in-plane boundary condition of the orbital 
anisotropy axis I at the A-B interface [Eq. ©]. The 
two calculated textures are the first with experimental 
basis, where fractionally quantized units appear in 3 He 
superfluids, since so far the predicted isolated i-quantum 
vortex ^3 has not been observed in 3 He-A. 

This collaboration was carried out under the EU-IHP 
ULTI-3, the ESF COSLAB, and ESF VORTEX pro- 
grams. 
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